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USE OF NUCLEOTIDE ANALOGS IN THE ANALYSIS OF 
OLIGONUCLEOTIDE MIXTURES AND IN HIGHLY MULTIPLEXED NUCLEIC 
ACID SEQUENCING 

Subject matter described herein was developed under NSF Grant 
No. Ger-9452651. The Government can have certain rights therein. 
RELATED APPLICATIONS 

This application is a divisional of allowed U.S. application Serial No. 
5 09/880,988, filed June 13, 2001, to Charles R. Cantor and Fouad A. 
Siddiqi, entitled "USE OF NUCLEOTIDE ANALOGS IN THE ANALYSIS OF 
OLIGONUCLEOTIDE MIXTURES AND IN HIGHLY MULTIPLEXED NUCLEIC 
ACID SEQUENCING." U.S. application Serial No. 09/880,988, claims 
the benefit of priority under 35 U.S.C. §1 19(e) to U.S. provisional 
10 application Serial No. 60/211,356, filed June 13, 2000, to Charles R. 

Cantor and Fouad A. Siddiqi, entitled USE OF NUCLEOTIDE ANALOGS IN 
THE ANALYSIS OF OLIGONUCLEOTIDE MIXTURES AND IN HIGHLY 
MULTIPLEXED NUCLEIC ACID SEQUENCING." Accordingly, the benefit 
of the priority date of U.S. provisional application Serial No. 60/21 1 ,356 
15 and U.S. application Serial No. 09/880,988 is claimed. The subject 

matter of each of U.S. provisional application Serial No. 60/21 1,356 and 
U.S. application Serial No. 09/880,988 is incorporated by reference in it 
entirety. 

FIELD OF THE INVENTION 
20 This invention relates to methods, particularly mass spectrometric 

methods, for the analysis and sequencing of nucleic acid molecules. 

DESCRIPTION OF THE BACKGROUND 

Since the recognition of nucleic acid as the carrier of the genetic 
code, a great deal of interest has centered around determining the 
25 sequence of that code in the many forms in which it occurs. Two studies 
made the process of nucleic acid sequencing, at least with DNA, a 
common and relatively rapid procedure practiced in most laboratories. 
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The first describes a process whereby terminally labeled DNA molecules 
are chemically cleaved at single base repetitions (A.M. Maxam and W. 
Gilbert, Proc. Natl. Acad. Sci. USA 74:560-64, 1977). Each base 
position in the nucleic acid sequence is then determined from the 
5 molecular weights of fragments produced by partial cleavage. Individual 
reactions were devised to cleave preferentially at guanine, at adenine, at 
cytosine and thymine, and at cytosine alone. When the products of these 
four reactions are resolved by molecular weight, using, for example, 
polyacrylamide gel electrophoresis, DNA sequences can be read from the 
10 pattern of fragments on the resolved gel. 

In another method DNA is sequenced using a variation of the plus- 
minus method (Sanger et at. (1977) Proc. Natl. Acad. Sci. USA 74:5463- 
67, 1977). This procedure takes advantage of the chain terminating 
ability of dideoxynucleoside triphosphates (ddNTPs) and the ability of 
15 DNA polymerase to incorporate ddNTPs with nearly equal fidelity as the 
natural substrate of DNA polymerase, deoxynucleoside triphosphates 
(dNTPs). Briefly, a primer, usually an oligonucleotide, and a template 
DNA are incubated in the presence of a useful concentration of all four 
dNTPs plus a limited amount of a single ddNTP. The DNA polymerase 
20 occasionally incorporates a dideoxynucleotide that terminates chain 
extension. Because the dideoxynucleotide has no 3'-hydroxyl, the 
initiation point for the polymerase enzyme is lost. Polymerization 
produces a mixture of fragments of varied sizes, all having identical 3' 
termini. Fractionation of the mixture by, for example, polyacrylamide gel 
25 electrophoresis, produces a pattern that indicates the presence and 

position of each base in the nucleic acid. Reactions with each of the four 
ddNTPs permits the nucleic acid sequence to be read from a resolved gel. 

These procedures are cumbersome and are limited to sequencing 
DNA. In addition, with conventional procedures, individual sequences are 
30 separated by, for example, electrophoresis using capillary or slab gels, 
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which stow. Mass spectrometry has been adapted and used for 
sequencing and detection of nucleic acid molecules (see, e.g., U.S. Patent 
Nos. (6,194,144; 6,225,450; 5,691,141; 5,547,835; 6,238,871; 
5,605,798; 6,043,031; 6,197,498; 6,235,478; 6,221,601; 6,221,605). 
In particular, Matrix-Assisted Laser Desorption/lonization (MALDI) and 
ElectroSpray Ionization (ESI), which allow intact ionization, detection and 
exact mass determination of large molecules, i.e. well exceeding 300 kDa 
in mass have been used for sequencing of nucleic acid molecules. 

A further refinement in mass spectrometric analysis of high 
molecular weight molecules was the development of time of flight mass 
spectrometry (TOF-MS) with matrix-assisted laser desorption ionization 
(MALDI). This process involves placing the sample into a matrix that 
contains molecules that assist in the desorption process by absorbing 
energy at the frequency used to desorb the sample. Time of flight 
analysis uses the travel time or flight time of the various ionic species as 
an accurate indicator of molecular mass. Due to its speed and high 
resolution, time-of-flight mass spectrometry is well-suited to the task of 
short-range, i.e., less than 30 base sequencing of nucleic acids. Since 
each of the four naturally occurring nucleotide bases dC, dT, dA and dG, 
also referred to herein as C, T, A and G, in DNA has a different molecular 
weight, 

Mc = 289.2 
Mt = 304.2 
Ma = 313.2 
Mg = 329.2, 

where Mc, Mt, Ma, Mg are average molecular weights in daltons of the 
nucleotide bases deoxycytidine, thymidine, deoxyadenosine, and 
deoxyguanosine, respectively, it is possible to read an entire sequence in 
a single mass spectrum. If a single spectrum is used to analyze the 
products of a conventional Sanger sequencing reaction, where chain 
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termination is achieved at every base position by the incorporation of 
dideoxynucleotides, a base sequence can be determined by calculation of 
the mass differences between adjacent peaks. In addition, the method 
can be used to determine the masses, lengths and base compositions of 
5 mixtures of oligonucleotides and to detect target oligonucleotides based 
upon molecular weight. 

MALDI-TOF mass spectrometry for sequencing DNA using mass 
modification (see, e.g., U.S. Patent Nos. 5,547,835, 6,194,144; 
6,225,450; 5,691,141 and 6,238,871) to increase mass resolution is 

10 available. The methods employ conventional Sanger sequencing reactions 
with each of the four dideoxynucleotides. In addition, for example for 
multiplexing, two of the four natural bases are replaced; dG is substituted 
with 7-deaza-dG and dA with 7-deaza-dA. 

A variety of techniques and combinations thereof have been 

15 directed to improving the level of accuracy in determining the nucleotide 
compositions of mixtures of oligonucleotides using mass spectrometry, 
and many of these methods employ nucleotide analogs. For example, 
Muddiman eta/. {Anal. Chem., 69(8): 1543-1549, 1997) discusses an 
algorithm for the unique definition of the base composition of PCR- 

20 amplified products, especially longer (> 100bp) oligonucleotides. The 
algorithm places a constraint on the otherwise large number of possible 
base compositions for long oligonucleotides by taking into account only 
those masses (measured by electrospray ionization mass spectrometry) 
that are consistent with that of their denatured complementary strands, 

25 assuming Watson-Crick base-pairing. In addition, the algorithm imposes 
the constraint of known primer compositions, since the primer sequences 
are known, and this constraint becomes especially significant with shorter 
PCR products whose mass of "unknown" sequence relative to that of the 
primer mass is small. Muddiman eta/, also discusses invoking additional 

30 measurements for defining the base composition with even greater 
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accuracy. These include the possibility of post-modifying the PCR 
product using e.g., dimethyl sulfate to selectively methylate every "G" in 
the PCR product, or using a modified base during PCR amplification, 
conducting mass measurements on the modified oligonucleotides, and 
5 comparing the mass measurements with those of the unmodified 
complementary strands. 

Chen eta/. (Anal. Chem., 71 (1 5): 31 1 8-31 25, 1999) reports a 
method that combines stable isotope 13 C/ 15 N labelling of PCR products 
with analysis of the mass shifts by MALDI-TOF mass spectrometry. The 

10 mass shift due to labelling of a single type of nucleotide (i.e, A, T, G or C) 
reveals the number of that type of nucleotide in a given fragment. While 
the method is useful in the measurement and comparison of nucleotide 
compositions of homologous sequences for sequence validation and in 
scoring polymorphisms, tedious repetitive sequencing reactions (using the 

15 four different labelled nucleotides) and mass spectrometric measurements 
are required. 

Hence there is a need in the art for methods that (i) unambiguously 
assign nucleotides in a sequence, and, (ii) resolve large numbers of 
oligonucleotides that have the same length, different base compositions, 

20 and nearly equal (i.e., less than or equal to about 1 dalton difference) 
molecular weights. Therefore it is an object herein to provide methods 
that solve such problems 
SUMMARY OF THE INVENTION 

Provided herein are methods for sequencing and detecting nucleic 

25 acids using techniques, such as mass spectrometry and gel 

electrophoresis, that are based upon molecular mass. The methods use 
deoxynucleotide analogs, modified nucleotide terminators and/or mass- 
labeled primers in one or more reactions for sequencing or detection 
protocols that involve primer extension, and analyze these products from 

30 more than one oligonucleotide on, for example, a single mass spectrum. 
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This provides a means for accurate detection and/or sequencing of a an 
oligonucleotide and is particularly advantageous for detecting or 
sequencing a plurality target nucleic acid molecules in a single reaction 
using any technique that distinguishes products based upon molecular 
5 weight. The methods herein are particularly adapted for mass 
spectrometric analyses. 

For example, a sequencing method provided herein uses 
deoxynucleotide analogs, modified nucleotide terminators and/or mass- 
labeled primers in one or more Sanger sequencing reactions, and analyzes 

10 these products from more than one oligonucleotide on a single mass 
spectrum. In particular, a plurality of primers can be used to 
simultaneously sequence a plurality of nucleic acid molecules or portions 
of the same molecule. By extending the primers with mass-matched 
nucleotides, the resulting products mass shifts that are periodically related 

15 to the size of the original primer. 

As a result, the sequence of any given oligonucleotide can be 
determined with a high level of accuracy, and also mixtures of a number 
of sequences can be multiplexed in a single mass spectrum. The limit on 
the number of oligonucleotides that can be sequenced simultaneously is 

20 governed by the base periodicity, the maximum mass shift, and the 

resolving power of analytical tool, such as the mass spectrometer. The 
base periodicity and maximum mass shift can be carefully engineered for 
optimal resolution and accuracy, depending on the number of sequences 
to be simultaneously analyzed, and the information desired; as many 

25 sequences as desired can be sequenced simultaneously especially in the 
detection and scoring of single nucleotide polymorphisms, insertions, 
deletions and other mutations. 

In another embodiment, a target nucleic acid molecule is 
sequenced using mass-matched nucleotides and chain terminating 

30 nucleotides. For example, a primer is annealed to a target nucleic acid, 
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the primer is extended in the presence of chain-terminating nucleotides 
and mass-matched nucleotides to produce extension products, the 
masses of the extension products follow a periodic distribution that is 
determined by the mass of the mass-matched nucleotides, and the 
5 sequence of the target nucleic acid is determined from the mass shift of 
each extension product from its corresponding periodic reference mass by 
virtue of incorporation of the chain terminator. The mass-matched 
nucleotides all have identical masses, and each chain terminating 
nucleotide has a distinct mass that differs from that of the other chain 

10 terminating nucleotides. This results in unique predetermined values of 
mass shift corresponding to each chain terminating nucleotide and based 
upon the original primer. 

This method is adaptable for any sequencing method or detection 
method that relies upon or includes chain extension. These methods 

15 include, but are not limited to, sequencing methods based upon Sanger 
sequencing, and detection methods, such as primer oligo base extension 
(PROBE) (see, e.g., U.S. application Serial No. 6,043,031; allowed U.S. 
application Serial No. 09/287,679; and 6,235,478), that rely include a 
step of chain extension. 

20 Also, contemplated are methods, such as haplotyping methods, in which 
two mutations in the same gene are detection are provided. A detector 
(primer) oligonucleotide is to the hybridized to the first mutation and the 
primer is extended with mass-matched nucleotides and appropriately 
selected chain terminator (s) to detect the second mutation. 

25 In other embodiments, a plurality of target nucleic acids can be 

multiplexed in a single reaction measurement by annealing each target 
nucleic acid to a primer of distinct molecular weight each primer is then 
extended with mass-matched nucleotides and chain terminators in 
formats that depend upon whether detection or sequencing is desired. 

30 These methods are particularly useful for methods of detection in which a 
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primer is hybridized to a plurality of target nucleic acid molecules, such as 
immobilized nucleic acid molecules, hybrids separated from unhybridized 
nucleic acids and the detectors detected. Such methods include PROBE, 
in which case the extension reaction is performed in the presence chain 
5 terminators and mass matched deoxynucleotides. 

The primers of distinct molecular weight can be selected to differ in 
molecular weight by a value that is greater than the maximum mass shift, 
i.e., the difference in molecular weight between the heaviest and the 
lightest nucleotide terminators in chain extension reactions. The 
10 difference in molecular weight between the primers for a plurality of 
target nucleic acids can be selected to be least 20 daltons greater than 
the maximum mass shift to account for the finite band width of the 
peaks. 

The number of molecules that can be multiplexed is governed by 

15 the periodicity, the maximum mass shift, and the resolving power of the 
sequence detection instrument. In some embodiments, about 7 to about 
25 or more molecules can be multiplexed. For scoring single nucleotide 
polymorphisms, only a single nucleotide terminator is required (depending 
on the base identity of the single nucleotide polymorphism). In this case, 

20 the maximum mass shift required is identically zero, so that larger 
numbers of molecules, greater than 25, 35, 50 and more, can be 
multiplexed, depending on the resolving power of the sequencing format, 
and for mass spectrometry the instrument. Depending on the amount of 
sequence information desired, one, two or three rather than four types of 

25 nucleotide terminators (corresponding to each of the four nucleic acid 
bases) can be used. 

In other embodiments, the mass shift is obtained using pair- 
matched nucleotides, i.e., the mass of each nucleotide base-pair is 
selected so that the masses of all pairs are identical. In one embodiment 

30 thereof, the following steps are performed: (i) the target nucleic acid is 
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copied or amplified by a method such as PCR in the presence of the pair- 
matched nucleotide set prior to the sequencing or detection reaction; 
(ii) the target nucleic acid is denatured, and a partially duplex hairpin 
primer is annealed and ligated to the single-stranded template; (iii) the 
5 primer is extended in the presence of chain terminating nucleotides and 
pair-matched nucleotides to produce extension products, where the 
masses of the extension products follow a periodic distribution that is 
determined by the mass of the pair-matched nucleotide set, and, (iv) the 
target nucleic acid is detected by virtue of its molecular weight or its 
10 sequence is determined from the mass shift of each extension product 
from its corresponding periodic reference mass. 

In another embodiment, the mass of each terminating base pair is 
unique and resolvable, so that the mass shifts corresponding to each 
terminating base pair are unique. The nucleotide terminators are 
15 optionally mass-matched or can be of distinct masses as long as distinct 
values of mass shift are obtained for each terminating base pair. 

In another embodiment, the extension products are treated to 
produce blunt-ended double-stranded extension products by methods 
known to those of skill in the art, such as the use of single-strand specific 
20 nucleases. In an aspect of this embodiment, a plurality of target nucleic 
acids can be multiplexed in a single reaction by annealing each target 
nucleic acid to a primer of distinct molecular weight. The primers can be 
selected to differ in molecular weight by a value that is greater than the 
maximum mass shift, i.e., the difference in molecular weight between the 
25 heaviest and the lightest nucleotide terminating base pairs. Since double 
stranded nucleic acid can be analyzed, the effective sequence read is 
halved relative to the embodiment employing mass-matched nucleotides, 
but the number of molecules that can be multiplexed is doubled, due to 
the increase in period (the value of the mass of a base pair, rather than a 
30 single mass-matched nucleotide). In exemplary embodiments, about 14 
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to about 50 sequences are multiplexed. In detection embodiments, about 
50 or more molecules can be simultaneously multiplexed since only a 
single terminating base pair is added in the extension reaction. 

In another embodiment, the chain termination reactions are carried 
5 out separately using a standard nucleotide terminator, pair-matched 

nucleotides, and mass-labeled primers, if modified nucleotide terminators 
which are either mass-matched or provide distinct values of mass shift for 
each terminating base pair are not available. The reactions are pooled 
prior to detection or sequence analysis. In one embodiment, the mass- 

10 labeled primers can have distinct values of molecular weight that give rise 
to unique values of mass shift or positional mass difference for each 
terminating base. 

In andother method provided herein, a population of nucleic acids 
having the same length but different base compositions can be resolved 

15 by synthesizing the nucleic acids in the presence of a nucleotide analog to 
produce synthesized nucleic acids having incorporated the nucleotide 
analog, where the nucleotide analog is selected to optimally separate the 
masses of the population of nucleic acids according to their individual 
base compositions. For example, the nucleotide analog or analogs are 

20 selected to separate the population of nucleic acids according to base 
composition by greater than 1 dalton. In another embodiment, the 
nucleotide analog or analogs are selected to separate the population of 
nucleic acids according to base composition by mass values of about 3 
daltons to about 8 daltons, depending on the choice of analog and on the 

25 resolving power of the detection instrument. In other embodiments, the 
nucleotide analog or analogs can be selected to restrict oligonucleotides 
having the same length to have the same mass, i.e., a peak separation of 
zero, regardless of differences in base composition, such as in detection 
methods, where it is desirable to separate populations of oligonucleotides 

30 according to their length. 
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Nucleic acid molecules that contain mass-matched nucleotides 
and/or pair-matched nucleotides are provided. 

Also provided are combinations for practicing the methods provided 
herein. For instance, in one embodiment, the combinations include a set 
5 of mass-matched deoxynucleotides. In another embodiment, the 
combinations a set of pair-matched nucleotides and a set of mass- 
matched chain terminating nucleotides. In another embodiment, the 
combination includes a set of pair-matched nucleotides and chain 
terminating nucleotides which form terminating base pairs of distinctly 
10 different molecular weight. In yet another embodiment, the combination 
includes a set of pair-matched nucleotides and mass-labeled primers. In 
other embodiments, mass-staggered primers can be added to as optional 
components. 

Kits containing the combinations with optional instructions and/or 
15 additional reagents are also provided. The kits contain the reagents as 

described herein and optionally any other reagents required to perform the 
reactions. Such reagents and compositions are packaged in standard 
packaging known to those of skill in the art. Additional vials, containers, 
pipets, syringes and other products for sequencing can also be included. 
20 Instructions for performing the reactions can be included. 

Also provided herein are methods for optimization of the analysis of 
base compositions of mixtures of oligonucleotides by mass spectrometry. 
A single spectrum can be used to resolve a very large number of 
oligonucleotides having the same length but different molecular weights 
25 by incorporating a nucleotide analog into the oligonucleotides in the 
mixture such that the peaks are no closer than a minimum value called 
peak separation. The peak separation can be tailored by careful selection 
of the nucleotide analog and of a mass spectrometer with the desired 
resolving power. 
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The methods herein permit unambiguous and accurate analysis of 
the sequences or molecular weights of large numbers of oligonucleotides 
in a single mass spectrum by combining the rapidity of mass spectrometry 
with the resolving power of nucleotide analogs which are carefully 
5 selected and incorporated into the oligonucleotide mixture according to 
the desired application. 

Other features and advantages will be apparent from the following 
detailed description and claims. 
BRIEF DESCRIPTION OF THE FIGURES 
10 Figure 1 shows that when a single spectrum is used to analyze the 

products of a conventional Sanger sequencing reaction, where chain 
termination is achieved at every base position by the incorporation of 
dideoxynucleotides, the base sequence can be determined by calculation 
of the mass differences between adjacent peaks (Figures 1a and 1b). 
15 Figure 2 shows implementation of forced mass modulation using 

mass-matched deoxynucleotides. Figure 2a is a simulated mass spectrum 
showing the products and molecular masses of a reaction carried out with 
a suitable polymerase in the presence of a mass-matched nucleotide set 
("dN") and the four standard dideoxynucleotide terminators. The base 
20 periodicity is the mass of dN, or 310 daltons. Figure 2b shows a target 
second sequence resolved on the same mass spectrum shown in Figure 
2a, using a primer heavier by 77 daltons. The peaks corresponding to the 
reaction products from the first target sequence can fall within the 
spectrum in Figure 2b, which can never intersect peaks from the second 
25 target sequence. This permits unambiguous resolution of both sequences 
each peak can be uniquely assigned to a nucleotide, a base position, and 
a target sequence. 

Figure 3 shows four different sequences resolved in a single 
spectrum using a set of mass-staggered primers that are separated in 
30 mass by integer multiples of 77 daltons (77, 154, and 231 daltons). 
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Figure 4 shows the general implementation of a forced mass 
modulation method using pair-matched nucleotides, for the analysis of 
sequencing reaction products as double-stranded structures. The steps in 
the reaction are as follows: a) a partially duplex hairpin primer with a 3' 
5 overhang and a 5' phosphate group is annealed and ligated to the single 
stranded target sequence; b) the resulting partially duplex structure is 
subjected to a sequencing reaction using the pair-matched nucleotide set 
described above along with the set of mass-matched terminators (ddM); 
c) products resulting from sequencing reaction b); and, d) the products c) 

10 from the sequencing reaction are exposed to a strict single strand-specific 
nuclease that results in the production of blunt-ended hairpin structures 
ready for analysis by mass spectrometry. 

Figure 5 shows the products and molecular masses of the nuclease 
digestion elucidated in Figure 4d, along with a simulated mass spectrum. 

15 Figure 6 shows three sequence variants (Figure 6a) that differ from 

each other only at a single base position sequenced by a conventional 
Sanger reaction. Figure 6b is a simulated mass spectrum of all reaction 
products shown in Figure 6a. Figure 6c is a graph representing the valid 
sequence permutations that can be elucidated from the mass spectrum 

20 shown in Figure 6b. Boxed values are fragment masses, solid arrows 
show valid sequence branches, dashed arrows represent spurious 
branches. In practice valid branches are indistinguishable from spurious 
ones. Figure 6d is a set of sequences consistent with the graph shown in 
Figure 6c. Spurious sequence reconstructions are shown in lowercase 

25 letters, valid ones in uppercase letters. 

Figure 7a shows the three related sequence variants from Figure 6 
sequenced by Forced Mass Modulation using a single primer and the 
mass-matched nucleotide set from Figure 2 with the standard dideoxy 
terminators. The positions of the differing bases are shown by solid 

30 arrows. Reaction products are shown along with their respective 
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molecular masses. Reaction products of variant #2 whose masses differ 
from those of variant #1 are marked with by (*). Reaction products of 
variant #3 whose masses differ from those of variant #1 are marked by 
(**). Figure 7b is a simulated mass spectrum of all reaction products 
shown in Figure 7a along with a sequence graph. The shaded regions 
represent the only valid mass ranges that can be assumed by the reaction 
products from Figure 7a. The base periodicity is 310 daltons. Figure 7c 
is a consensus sequence derived from the data shown in Figure 7b. 
Figure 7d is an expansion of the consensus sequence shown in Figure 7c. 
Spurious reconstructions are shown in lowercase letters, valid ones in 
uppercase letters. Note that there is only a single spurious 
reconstruction, as opposed to the eleven errant sequences reconstructed 
from the Sanger reaction described in Figure 6. 

Figure 8 shows the base composition density distributions for the 
total set of possible 7-base oligonucleotides using three different 
nucleotide sets. Note that for the set of naturally occurring bases, nearly 
every base composition has its own distinct mass value, but most of 
these mass values are spaced only one dalton from each other. 
Increasing the peak separation also markedly increases the average 
number of base compositions per observed mass, particularly for those 
masses in the center of the range. 
DETAILED DESCRIPTION OF THE INVENTION 
Definitions 

Unless defined otherwise, all technical and scientific terms used 
herein have the same meaning as is commonly understood by one of skill 
in the art to which this invention belongs. All patents, patent 
applications, Genbank and other sequence repository sequences, and 
publications referred to herein are incorporated by reference. 

As used herein, a biopolymer includes, but is not limited to, nucleic 
acid, proteins, polysaccharides, lipids and other macromolecules. Nucleic 

-14- 



25491 -2408B 



acids include DNA, RNA, and fragments thereof. Nucleic acids may be 
derived from genomic DNA, RNA, mitochondrial nucleic acid, chloroplast 
nucleic acid and other organelles with separate genetic material. 

As used herein "nucleic acid" refers to polynucleotides such as 
5 deoxyribonucleic acid (DNA) and ribonucleic acid (RNA). The term should 
also be understood to include, as equivalents, derivatives, variants and 
analogs of either RNA or DNA made from nucleotide analogs, single 
(sense or antisense) and double-stranded polynucleotides. 
Deoxyribonucleotides include deoxyadenosine, deoxycytidine, 

10 deoxyguanosine and deoxythymidine. For RNA, the uracil base is uridine. 
As used herein, "forced mass modulation" refers to methods 
provided herein that use deoxynucleotide analogs, modified nucleotide 
terminators, mass-labeled primers, mass-staggered primers and other such 
nucleotides, nucleic acids and analogs thereof, to unambiguously assign 

15 peak positions of mass fragments of oligonucleotides according to their 
base position, base identity, and target sequence from which the 
fragments arose. The method is used to sequence, detect or identify 
single oligonucleotide or plurality thereof. Hence the method is used, for 
example for muliplex sequencing and detection of nucleic acid molecules 

20 among mixtures thereof. 

As used herein, "nucleotides" include, but are not limited to, the 
naturally occurring nucleoside mono-, di-, and triphosphates: 
deoxyadenosine mono-, di- and triphosphate; deoxyguanosine mono-, di- 
and triphosphate; deoxythymidine mono-, di- and triphosphate; and 

25 deoxycytidine mono-, di- and triphosphate (referred to herein as dA, dG, 
dT and dC or A, G, T and C, respectively). Nucleotides also include, but 
are not limited to, modified nucleotides and nucleotide analogs such as 
deazapurine nucleotides, e.g., 7-deaza-deoxyguanosine (7-deaza-dG) and 
7-deaza-deoxyadenosine (7-deaza-dA) mono-, di- and triphosphates, 

30 deutero-deoxythymidine (deutero-dT) mon-, di- and triphosphates, 
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methylated nucleotides e.g., 5-methyldeoxycytidine triphosphate, 13 C/ 15 N 
labelled nucleotides and deoxyinosine mono-, di- and triphosphate. For 
those skilled in the art, it will be clear that modified nucleotides and 
nucleotide analogs can be obtained using a variety of combinations of 
functionality and attachment positions. 

As used herein, a complete set of chain-elongating nucleotides 
refers to the four different nucleotides or analogs thereof that hybridize to 
each of the four different bases comprising the nucleic acid template. 

As used herein, the term "mass-matched nucleotides" refers to a 
set of nucleotide analogs wherein each analog is of identical mass to each 
of the other analogs. For example, analogs of dA, dG, dC and dT can 
form a mass-matched nucleotide set, when each analog is selected to 
have the same molecular weight as the others in the set. Mass-matched 
nucleotide sets can be identified by selecting chemically modified 
derivatives of natural bases or by the use of a universal base analog such 
as deoxyinosine or 5-nitroindole and 3-nitropyrrole (5-nitroindole and 3- 
intropyrrole can be in the dideoxy form) which can form base pairs with 
more than one of the natural bases. Others include, 3-methyl 7-propynyl 
isocarbostyril, 5-methyl isocarbostyril, and 3-methyl isocarbostyril. As a 
result, oligonucleotides that contain such bases differ in molecular weight 
only as a function of length thereof. Furthermore, incorporation of a single 
nucleotide(s) that is (are) not in the set renders such the 
oligonucleotide(s) readily identifiable by mass, particularly by 
spectrometric analysis. 

As used herein, the term "pair-matched nucleotides" refers to a 
nucleotide set in which the nucleotide analogs are selected such that the 
total mass each base pair is identical. For example, replacing dG with the 
nucleotide analog 7-deaza-dG forces the mass of each base pair, i.e., (dA 
+ dT) and (dC + 7-deaza-dG) to be identical. Exemplary pair-matched 
nucleotides, include, but are not limited to, 7-deaza-dA + 
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phosphorothioate-dT ((312.2 + 320.2) = 632.4 Da) and 5-methyl-dC + 
dG ((303.2 + 329.2) = 632.4 Da); phosphorothioate-7-deaza-dA + dU 
((328.2 + 290.2) = 618.4 Da) and dC + dG = {(289.2 + 329.2) = 
618.4 Da), and other such pairs that may be readily selected. Another 
5 exemplary set of mass-matched nucleotides with a molecular mass of 
328.2: 7-deaza-dG, phosphorothioate-7-deaza-dA, 5-propynyl-dU and 5- 
cyanomethyl-2'-deoxycytidine. 

As used herein, the term "nucleotide terminator" or "chain 
terminating nucleotide" refers to a nucleotide analog that terminates 

10 nucleic acid polymer (chain) extension during procedures wherein a DNA 
template is being sequenced or replicated. The standard chain 
terminating nucleotides, i.e., nucleotide terminators include 2', 3'- 
dideoxynucleotides (ddATP, ddGTP, ddCTP and ddTTP, also referred to 
herein as dideoxynucleotide terminators). As used herein, 

15 dideoxynucleotide terminators also include analogs of the standard 

dideoxynucleotide terminators, e.g., 5-bromo-dideoxyuridine, 5-methyl- 
dideoxycytidine and dideoxyinosine are analogs of ddTTP, ddCTP and 
ddGTP, respectively. 

As used herein, "mass-matched terminators" refers to a set of 

20 nucleotide terminators that are selected such that each analog of ddA, 

ddG, ddC and ddT making up the mass-matched set has exactly the same 
molecular weight. Mass-matched terminator sets can be constructed by 
selecting chemically modified derivatives of standard dideoxynucleotides 
or by the use of a universal dideoxynucleotide analog that form base pairs 

25 with more than one of the natural bases. Exemplary mass-matched 
nucleotides include, but are not limited to, 3-methyl 7-propynyl 
isocarbostyril, 5-methyl isocarbostyril and 3-methyl isocarbostyril. As 
used herein, the terms "oligonucleotide" or "nucleic acid" refer to single- 
stranded and/or double-stranded polynucleotides such as deoxyribonucleic 

30 acid (DNA), and ribonucleic acid (RNA) as well as derivatives of either 
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RNA or DNA into which nucleotide or dideoxynucleotide analogs have 
been incorporated. Also included in the term "nucleic acid" are analogs 
of nucleic acids such as peptide nucleic acid (PNA), phosphorothioate 
DNA, and other such analogs and derivatives. 
5 As used herein, "nucleotide composition" or "base composition" 

refers to the numerical ratio of the four nucleotide bases relative to each 
other in an oligonucleotide. 

As used herein, a target nucleic acid refers to any nucleic acid of 
interest in a sample. It can contain one or more nucleotides. A target 

10 nucleotide sequence refers to a particular sequence of nucleotides in a 

target nucleic acid molecule. Detection or identification of such sequence 
results in detection of the target and can indicate the presence or absence 
of a particular mutation or polymorphism. 

As used herein, "partially duplex hairpin" refers to a partially self- 

15 complementary oligonucleotide, which forms intramolecular base-pairs 

within its self-complementary region, leaving a "loop" of bases at one end 
of the molecule and a single-stranded "overhang" region at the other end. 
Thus, the oligonucleotide assumes a hairpin-like motif. "Blunt-ended 
hairpin structures", as referred to herein, are similar to the partially duplex 

20 hairpin structures with the exception that they do not have a single- 
stranded "overhang" region. 

As used herein, "base periodicity" or "period" (Pbase) refers to the 
quasi-periodic distribution of the molecular weights of products obtained 
using Forced Mass Modulation. The base periodicity results from either 

25 the mass of the mass-matched deoxynucleotide set, or the mass of the 
pair-matched deoxynucleotide set, or from the modified chain terminators 
depending on the embodiment implemented. The base sequence or 
nucleic acid molecule identity is encoded in the pattern (or detectable 
therein) in which the observed mass distribution deviates from absolute 

30 regular periodicity. 
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As used herein, the "periodic reference mass" at base position "n" 
in any given oligonucleotide molecule, MpR[n], is defined as the sum of: (i) 
the mass of the primer (Mprimer) used to sequence the DNA template using 
Forced Mass Modulation, (ii) the mass of the lightest nucleotide 
5 terminator (Mnght), and, (iii) (n-1) multiple of the base periodicity Pbase. 

As used herein, the "positional mass difference" or "mass shift" at 
base position "n" in any given oligonucleotide molecule, Mdiff[n], is 
defined as the distance in daltons between the observed peak, Mobs[n], 
and the nth periodic reference mass. 

10 As used herein, the "maximum mass shift" Smax is the maximum 

possible value of the positional mass difference, depending on the choice 
of mass-matched nucleotides and nucleotide terminators used in the 
implementation of Forced Mass Modulation. Accordingly, the maximum 
mass shift can be modulated by the choice of mass-matched nucleotides 

15 and nucleotide terminators. 

As used herein, a "primer" refers to an oligonucleotide that is 
suitable for hybridizing, chain extension, amplification and sequencing. 
Similarly, a probe is a primer used for hybridization. The primer refers to 
a nucleic acid that is of low enough mass, typically about between about 

20 5 and 200 nucleotides, generally about 70 nucleotides or less than 70, 
and of sufficient size to be conveniently used in the methods of 
amplification and methods of detection and sequencing provided herein. 
These primers include, but are not limited to, primers for detection and 
sequencing of nucleic acids, which require a sufficient number nucleotides 

25 to form a stable duplex, typically about 6-30 nucleotides, about 10-25 

nucleotides and/or about 12-20 nucleotides. Thus, for purposes herein, a 
primer is a sequence of nucleotides of any suitable length, typically 
containing about 6-70 nucleotides, 12-70 nucleotides or greater than 
about 14 to an upper limit of about 70 nucleotides, depending upon 

30 sequence and application of the primer. 
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As used herein, the term "mass-labeled primers" refers to a set of 
primers that differ in mass by values that provide distinct and resolvable 
positional mass differences for each of the four termination reactions in 
an embodiment of Forced Mass Modulation. In this particular 
5 embodiment of Forced Mass Modulation, each of the termination 

reactions for a given oligonucleotide is carried out separately using each 
of the mass-labeled primers, and the reaction products are combined prior 
to obtaining a mass spectrum. 

As used herein, the term "mass-staggered primers" refers to the 
10 mass difference ("staggering" of the masses) between the primers used 
in multiplexed sequencing using Forced Mass Modulation. For resolution 
of multiple sequences using this method, the differences between the 
masses of the primers should at least be equal to the maximum mass 
shift, and is generally greater than the maximum mass shift by at least 20 
15 daltons to account for the finite width of each observed peak. 

As used herein, reference to mass spectrometry encompasses any 
suitable mass spectrometric format known to those of skill in the art. 
Such formats include, but are not limited to, Matrix-Assisted Laser 
Desorption/lonization, Time-of-Flight (MALDI-TOF), Electrospray (ES), IR- 
20 MALDI (see, e.g., published International PCT application No. 99/5731 8 
and U.S. Patent No. 5,118,937), Ion Cyclotron Resonance (ICR), Fourier 
Transform and combinations thereof. MALDI, particular UV and IR, are 
among the preferred formats. 

As used herein, mass spectrum refers to the presentation of data 
25 obtained from analyzing a biopolymer or fragment thereof by mass 
spectrometry either graphically or encoded numerically. 

As used herein, pattern with reference to a mass spectrum or mass 
spectrometric analyses, refers to a characteristic distribution and number 
of signals (such peaks or digital representations thereof). 
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As used herein, signal in the context of a mass spectrum and 
analysis thereof refers to the output data, which the number or relative 
number of molecules having a particular mass. Signals include "peaks" 
and digital representations thereof. 
5 As used herein, "mass spectrum division multiplexing" is an 

embodiment of Forced Mass Modulation in which unambiguous resolution 
of multiple sequences in a single spectrum is possible by judicious 
selection of mass staggered primers. 

As used herein, "analysis" refers to the determination of certain 

10 properties of a single oligonucleotide, or of mixtures of oligonucleotides. 
These properties include, but are not limited to, the nucleotide 
composition and complete sequence of an oligonucleotide or of mixtures 
of oligonucleotides, the existence of single nucleotide polymorphisms 
between more than one oligonucleotide, the masses and the lengths of 

15 oligonucleotides and the presence of a molecule or sequence within 
molecule in a sample. 

As used herein, "multiplexing" refers to the simultaneous 
determination of more than one oligonucleotide molecule, or the 
simultaneous analysis of more than one oligonucleotide, in a single mass 

20 spectrometry or other sequence measurement, i.e., a single mass 
spectrum or other method of reading sequence. 

As used herein, "polymorphisms" refer to variants of a gene or an 
oligonucleotide molecule that differ at more than one base position. In 
"single nucleotide polymorphisms", the variants differ at only a single 

25 base position. 

As used herein, amplifying refers to means for increasing the 
amount of a bipolymer, especially nucleic acids. Based on the 5' and 3' 
primers that are chosen, amplification also serves to restrict and define 
the region of the genome which is subject to analysis. Amplification can 

30 be by any means known to those skilled in the art, including use of the 
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polymerase chain reaction (PCR) etc. Amplification, e.g., PCR must be 
done quantitatively when the frequency of polymorphism is required to be 
determined. 

As used herein, "polymorphism" refers to the coexistence of more 
5 than one form of a gene or portion thereof. A portion of a gene of which 
there are at least two different forms, i.e., two different nucleotide 
sequences, is referred to as a "polymorphic region of a gene". A 
polymorphic region can be a single nucleotide, the identity of which 
differs in different alleles. A polymorphic region can also be several 

10 nucleotides in length. Thus, a polymorphism, e.g. genetic variation, 

refers to a variation in the sequence of a gene in the genome amongst a 
population, such as allelic variations and other variations that arise or are 
observed. Thus, a polymorphism refers to the occurrence of two or more 
genetically determined alternative sequences or alleles in a population. 

15 These differences can occur in coding and non-coding portions of the 

genome, and can be manifested or detected as differences in nucleic acid 
sequences, gene expression, including, for example transcription, 
processing, translation, transport, protein processing, trafficking, DNA 
synthesis, expressed proteins, other gene products or products of 

20 biochemical pathways or in post-translational modifications and any other 
differences manifested amongst members of a population. A single 
nucleotide polymorphism (SNP) refers to a polymorphism that arises as 
the result of a single base change, such as an insertion, deletion or 
change in a base. 

25 A polymorphic marker or site is the locus at which divergence 

occurs. Such site may be as small as one base pair (an SNP). 
Polymorphic markers include, but are not limited to, restriction fragment 
length polymorphisms, variable number of tandem repeats (VNTR's), 
hypervariable regions, minisatellites, dinucleotide repeats, trinucleotide 

30 repeats, tetranucleotide repeats and other repeating patterns, simple 
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sequence repeats and insertional elements, such as Alu. Polymorphic 
forms also are manifested as different mendelian alleles for a gene. 
Polymorphisms may be observed by differences in proteins, protein 
modifications, RNA expression modification, DNA and RNA methylation, 
5 regulatory factors that alter gene expression and DNA replication, and any 
other manifestation of alterations in genomic nucleic acid or organelle 
nucleic acids. 

As used herein, "polymorphic gene" refers to a gene having at least 
one polymorphic region. 

10 As used herein, "allele", which is used interchangeably herein with 

"allelic variant" refers to alternative forms of a gene or portions thereof. 
Alleles occupy the same locus or position on homologous chromosomes. 
When a subject has two identical alleles of a gene, the subject is said to 
be homozygous for the gene or allele. When a subject has two different 

15 alleles of a gene, the subject is said to be heterozygous for the gene. 

Alleles of a specific gene can differ from each other in a single nucleotide, 
or several nucleotides, and can include substitutions, deletions, and 
insertions of nucleotides. An allele of a gene can also be a form of a gene 
containing a mutation. 

20 As used herein, "predominant allele" refers to an allele that is 

represented in the greatest frequency for a given population. The allele or 
alleles that are present in lesser frequency are referred to as allelic 
variants. 

As used herein, a subject, includes, but is not limited to, animals, 
25 plants, bacteria, viruses, parasites and any other organism or entity that 
has nucleic acid. Among subjects are mammals, preferably, although not 
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necessarily, humans. A patient refers to a subject afflicted with a disease 
or disorder. 

As used herein, a phenotype refers to a set of parameters that 
includes any distinguishable trait of an organism. A phenotype can be 
5 physical traits and can be, in instances in which the subject is an animal, 
a mental trait, such as emotional traits. 

As used herein, "resolving power" of a mass spectrometer is the 
ion separation power of the instrument, i.e., it is a measure of the ability 
of the mass spectrometer to separate peaks representing different 
10 masses. The resolving power R is defined as m/Am, where m is the ion 
mass and Am is the difference in mass between two resolvable peaks in a 
mass spectrum. 

As used herein, "assignment" refers to a determination that the 
position of a nucleic acid fragment indicates a particular molecular weight 
15 and a particular terminal nucleotide. 

As used herein, "a" refers to one or more. 

As used herein, "plurality" refers to two or more, up to an amount 
that is governed by the base periodicity, the maximum mass shift, and the 
resolving power of the mass spectrometer. 

20 As used herein, an array refers to a pattern produced by three or 

more items, such as three or more loci on a solid support. 

As used herein, "distinct" refers to a unique value of molecular 
weight, mass shift or period that is different from every other value of 
molecular weight, mass shift or period in the measurement. 

25 As used herein, "unambiguous" refers to the unique assignment of 

a particular oligonucleotide fragment according to the identity of its 
terminal base position and, in the event that a number of molecules are 
multiplexed, that the peak representing an oligonucleotide fragment can 
also be uniquely assigned to a particular molecule. 
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As used herein, the symbols Mc, Mt, Ma and Mg are average 
molecular weights in daltons of the nucleotides deoxycytidine, thymidine, 
deoxyadenosine and deoxyguanosine, respectively, or of analogs thereof. 
Mavg, the average molecular weight of any given oligonucleotide is a 
5 function of the average molecular weights of each of the nucleotides 
comprising the oligonucleotide, the numbers c, t, a and g of each 
nucleotide present in the oligonucleotide, the length of the oligonucleotide 
n' that is the sum of c, t, a and g, and the constant k that represents the 
mass of any other chemical groups on the molecule, such as terminal 
10 phosphates. 

As used herein, Ntotal is the total number of possible base 
compositions for an oligonucleotide of length n'. 

As used herein, "peak separation" or "minimum peak separation" 
S refers to the minimum value of the distance between consecutive peaks 
15 in a mass spectrum that resolves a large number of oligonucleotides 
having the same lengths but different molecular weights, i.e., different 
base compositions. The peak separation, which can be tailored by careful 
selection of the nucleotide analogs incorporated into the oligonucleotide 
and by a mass spectrometer of desired resolving power, is usually a 
20 positive integer greater than one, and typically a positive integer greater 
than or equal to 3. For two oligonucleotides having the same length n' 
but different base compositions, their molecular weights will either 
correspond to the same peak if the molecular weights are identical, or to 
two peaks separated at least by a value equal to the peak separation. 
25 As used herein, L is the maximum number of allowed 

oligonucleotide masses for a given nucleotide set. It is directly 
proportional to the oligonucleotide length n' and the mass difference 
between the heaviest and lightest nucleotides in the set, and is inversely 
proportional to the peak separation. 
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As used herein, D refers to the average density of different base 
compositions per allowed mass value, given the set of all possible base 
compositions of an oligonucleotide of length n'. 

As used herein, Mheaw refers to the mass of the heaviest nucleotide, 
5 nucleotide terminator or terminating base pair in daltons, depending on 
the specific embodiment of Forced Mass Modulation being described. 

As used herein, Mught refers to the mass of the lightest nucleotide, 
nucleotide terminator or terminating base pair in daltons, depending on 
the specific embodiment of Forced Mass Modulation being described. 
10 As used herein, Mpdmer is the mass of the primer in daltons. 

As used herein, Mobs[n] is the observed mass of the sequencing 
reaction at the nth base position. 

As used herein, Mterm[n] refers to the mass in daltons of the nth 
terminating nucleotide. 
15 As used herein, L' is the theoretical upper limit on the number of 

sequences that be multiplexed in a single mass spectrum. L' is directly 
proportional to the base periodicity Pbase, and is inversely proportional to 
the maximum mass shift Smax. 

As used herein, Mdupiex is the mass in daltons of the fully duplex 
20 hairpin primer in the implementation of Forced Mass Modulation using 
pair-matched nucleotides. 

As used herein, MddM is the mass in daltons of a dideoxy terminator 
that belongs to a set of mass-matched terminators. 

As used herein, Mtarg[n] is the mass of the nth nucleotide past the 
25 priming site in the 3' to 5' direction in the target sequence, i.e., the 
oligonucleotide whose sequence is being determined. 

As used herein, "specifically hybridizes" refers to hybridization of a 
probe or primer only to a target sequence preferentially to a non-target 
sequence. Those of skill in the art are familiar with parameters that affect 
30 hybridization; such as temperature, probe or primer length and 
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composition, buffer composition and salt concentration and can readily 
adjust these parameters to achieve specific hybridization of a nucleic acid 
to a target sequence. 

As used herein, a biological sample refers to a sample of material 
obtained from or derived from biological material, such as, but are not 
limited to, body fluids, such blood, urine, cerebral spinal fluid and synovial 
fluid, tissues and organs. Derived from means that sample can be 
processed, such as by purification or isolation and/or amplification of 
nucleic acid molecules. 

As used herein, a composition refers to any mixture. It may be a 
solution, a suspension, liquid, powder, a paste, aqueous, non-aqueous or 
any combination thereof. 

As used herein, a combination refers to any association between 
two or among more items. 

As used herein, "kit" refers to a package that contains a 
combination and optionally instructions and/or reagents and apparatus for 
use with the combination. 

Forced Mass Modulation for analysis of nucleic acid molecules 
Time of flight analysis and drawbacks thereof 

While time-of-flight mass spectrometry offers a number of 
advantages over conventional techniques such as gel electrophoresis, the 
peculiar relationship between the masses of the bases in DNA 
complicates the analysis of complex mixtures of oligonucleotides by mass 
spectrometry. For a given oligonucleotide, the average molecular weight, 
Mavg, is given by the following equation: 

(i) Mavg = k + cMc + tlVh- + alvU + glvlc 
where Mc, Mt, Ma, Mg are the average molecular weights of each of the 
four nucleotide bases (cytosine, thymine, adenine, guanine) and c, t, a, g 
represent the number of each base present in the oligonucleotide. The 
term k is a constant representing the mass of any other chemical groups 
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on the molecule, such as terminal phosphates. Rearranging equation (i) to 
give the average molecular weight as a function of the length of the 
oligonucleotide in bases yields 

(ii) Mavg = k + n'Mc + t(MT - Mc) + a(M A - Mc) + g(Mc - Mc) 
where n', the oligonucleotide length, is defined as 

n'=c+t + a + g 
Substituting the masses of the naturally occurring bases in DNA (to one- 
tenth dalton): 

Mc = 289.2 
Mt = 304.2 
Ma = 313.2 
Mg = 329.2 

into equation (ii) yields 

(iii) Mavg = k + 289. 2n' + t(304.2 - 289.2) + a(313.2 - 289.2) 
+ g(329.2 - 289.2), 

which can be simplified to 

(iv) Mavg = k + 289. 2n' + 1 5t + 24a + 40g 

Close inspection of equation (iv) reveals that it is almost always possible 
to find two oligonucleotides of the same length but of different base 
composition whose average masses differ by only one dalton. For 
example, all 7-mers having a base composition of A2C2G2T have an 
average molecular weight of (2167.4 + k), while all 7-mers with the base 
composition A3CGT2 have an average molecular weight of (2166.4 + k). 
Since the following relation 

(Mc + Mg) = (Mt + Ma) + 1 
is always true for the naturally occurring bases in DNA, simply replacing 
one C and one G in an oligonucleotide with one A and one T will produce 
a new oligonucleotide exactly one dalton lighter. Many other "single- 
dalton difference" relations, such as 

4Ma = (Mc + Mt + 2Mg) + 1 
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can readily be found for the naturally occurring bases. 

Thus, the possibility always exists that two or more 

oligonucleotides of same length and different molecular weight (and, 

therefore, different base composition) will be too close in mass to be 

resolved by a time-of-flight instrument. Two oligonucleotides of same 

length but different molecular weight differ in base composition unless 

they are each composed of different nucleotide analogs, whereas two 

oligonucleotides of same length and same molecular weight can have 

either the same or different base compositions. This problem becomes 

increasingly severe with increasing oligonucleotide size, since the total 

number of possible base compositions, Ntotal, scales as a cubic function 

of the oligonucleotide length n', in bases: 

(v) Ntotal = (n' + 1)(n' + 2)(n' + 3) 
6 

The use of time-of-flight mass spectrometry in sequencing 
applications also poses several potential problems. The great drawback 
of sequencing by the Sanger method is that the molecular weights of the 
Sanger reaction products can appear virtually anywhere on the mass axis 
depending on the particular sequence being examined. As a result, the 
absolute mass of any single Sanger fragment has to 
be measured with sufficient accuracy to calculate its distance from the 
masses from the fragments above and below it. Thus, determination of 
the identity of a single base depends on the accuracy of two separate 
mass measurements. Any error in a determination mass of a single 
fragment affects the accuracy of two bases in the sequence. 

For longer sequences (30-50 bases), it may not be possible to 
determine the mass difference between adjacent peaks with sufficient 
accuracy to unambiguously determine base identity. This is particularly a 
problem for the nucleotides A and T, which differ in mass only by nine 
daltons. The problem is addressed by resolving each of the four 



-29- 



25491 -2408B 



termination reactions in a separate mass spectrum. In this case each 
peak functions essentially as binary signal indicating the presence of a 
base at a particular position, much as in conventional electrophoretic 
sequencing. Using separate spectra, however, increases read accuracy 
5 but at the expense of increasing the number of required mass 
measurements by a factor of four. 

It is possible to resolve two target sequences by the Sanger 
method in a single mass spectrum, provided that all products of the 
sequencing reactions have unique and resolvable masses, and multiplex 

10 methods using mass modified bases have been developed. But, where 
two or more reaction products have the same mass, then unambiguous 
reconstruction of the two target sequences is not possible (see, e.g., 
Figures 1c-e). In addition, there is no way to determine a priori which 
observed masses belong to a particular sequence. In practice, this means 

15 that multiplexed Sanger sequencing by mass spectrometry can be 

difficult. The methods provided herein resolve these problem and provide 
a way to determine a priori which masses are associated with extension 
of a particular primer. 

Forced Mass Modulation 

20 As noted above. Forced Mass Modulation refers to methods 

provided herein that permit unambiguous assignment of peak positions (or 
masses) to mass fragments of oligonucleotides according to their base 
position, base identity, and target sequence from which the fragments 
arose. The methods use deoxynucleotide analogs, modified nucleotide 

25 terminators, mass-labeled primers, mass-staggered primers and other such 
nucleotides, nucleic acids and analogs thereof to provide a means for 
deconvoluting complex mass spectra or output from other mass 
determining techniques. These methods permit deconvolution of highly 
multiplexed nucleic acid reaction mixtures for sequencing methods and 

30 detection methods that include a step of primer extension. In practicing 
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these methods, primers are extended using mass-matched nucleotides 
and chain terminators (or in some embodiments mass where it is only 
necessary to detect incorporation (or the absence of incorporated) mass- 
matched terminators and optionally mass-matched chain extending 
5 nucleotides). Because the sequence and/or molecular mass of a primer is 
known, and the extended nucleotides have the same molecular mass, a 
periodicity in molecular mass that is a function of molecular weight of the 
selected mass matched nucleotide(s) results. 

As described in more detail below, for sequencing reactions using 

10 chain terminators, the deviation from the periodicity results from 

incorporation of a chain terminator. The deviation is a function of the 
particular terminator incorporated. For detection methods, incorporation 
of a terminator will indicate the presence of a mutation (if the terminator 
is selected to pair with the first mutated nucleotide). Any shift from 

15 periodicity will indicate the presence of the mutation. These methods, 
thus provide a simple, reliable way to detect the presence of a mutations 
or target nucleotide(s) in a sequence and to sequence nucleic acids. 
Forced mass modulation can be used with any method, such as mass 
spectrometry and gel electrophoresis, that relies on molecular weight as 

20 an output. Mass spectrometry is exemplified herein. 

The methods, designated Forced Mass Modulation methods, 
provided herein, are implemented by suitable selection of nucleotides 
and/or chain terminators, such as by the use of deoxynucleotide analogs, 
modified nucleotide terminators and mass-labeled primers in one or more 

25 reactions. Forced Mass Modulation can be used to simultaneously 
sequence or detect large numbers of oligonucleotides, such as, for 
example, twenty-five or more, with a high degree of resolution and 
accuracy. It can also be used to simplify the analysis of closely related 
sequence variants, as is required in the detection and scoring of 

30 nucleotide polymorphisms, including single nucleotide polymorphism 
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(SNPs) and for other genotypical analyses. Forced Mass Modulation 
greatly improves the use of mass spectrometry for nucleic acid analyses. 
Nearly every application relies on mass measurements that can benefit in 
increased accuracy and in a reduction of the number of required spectra. 
5 Another advantage of Forced Mass Modulation is the number of different 
ways in which it can be implemented, allowing it to be tailored to 
particular experimental or instrumental limitations. 

For example, compared to the conventional Sanger methods, 
Forced Mass Modulation, provides increased accuracy, simplified 

10 interpretation of mass data, and the ability to use a single mass spectrum 
for the unambiguous resolution of several distinct nucleic acid molecules. 
For mass spectrometry applications, the methods provide unambiguous 
assignment of peak positions of mass fragments of oligonucleotides 
according to their base position, base identity, and target sequence from 

15 which the fragments arose. Thus, the methods herein are 

advantageously used for multiplexing, in which a plurality of reactions are 
run in a single reaction (single pot). Forced Mass Modulation, exemplified 
with reference to sequencing methods, such as PROBE, can also be 
adapted to detection methods in which a primer is extended. 

20 In Forced Mass Modulation in which a primer is extended with 

mass-matched nucleotides, for examples, the molecular weights of 
extended nucleic acid chains, such as sequencing reaction products, are 
constrained since all extension products from the same primer will have a 
molecular weight that differs either by the length of the extension and the 

25 chain terminator. As a result, the extension products assume a 

quasi-periodic distribution on the mass axis with a predetermined base 
periodicity. For sequencing, the base sequence itself is encoded in the 
pattern in which the observed mass distribution deviates from absolute 
regular periodicity. Since the base periodicity will always be known a 

30 priori, since the primer is known, each peak in the observed mass 
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spectrum can be matched unambiguously to a unique nucleotide position 
in the target sequence. The initiating primers fix each set of nested 
fragments or extended products, and the use of mass-matched 
nucleotides creates the periodicity. 
5 As demonstrated by the Examples below, the method is 

advantageous for numerous applications including sequencing and a 
variety of detection methods, including primer oligo base extension 
(PROBE) {see, e.g., U.S. Patent No. 6,043,031; allowed U.S. application 
Serial No. 09/287,679, now U.S. Patent No. 6,258,538; and U.S. Patent 

10 No. 6,235,478) that use mass spectrometry to distinguish between 

extended primers. If the base compositions of the target oligonucleotides 
are known a priori then it is possible to select a nucleotide set that 
produces oligonucleotide masses that are distinct and resolvable for any 
particular instrument or application. 

15 Conversely, it is also possible to select a nucleotide set that 

restricts specific oligonucleotides to have the same mass, regardless of a 
change in base composition. The strategy of restricting specific 
oligonucleotides to have the same mass can be used to separate more 
than one oligonucleotide population of different lengths by restricting all 

20 oligonucleotides of a particular length to the same molecular weight, 
irrespective of differences in base composition. 

The oligonucleotide analysis or sequencing in methods provided 
herein can be accomplished by one of several methods employed in the 
art for the synthesis, resolution and/or detection of nucleic acids. 

25 Depending on the embodiment implemented, modified nucleotides can be 
incorporated into the oligonucleotides by chemical {Oligonucleotides and 
Analogues: A Practical Approach, F. Eckstein, ed., IRL Press Oxford, 
1991) or enzymatic (F. Sanger et al., Proc. Natl. Acad. Set. USA 
74:5463-67, 1977) synthesis. Extension products or truncated products 

30 of the oligonucleotides to be sequenced can be obtained using chemical 
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(A.M. Maxam and W. Gilbert, Proc. Natl. Acad. Sci. USA 74:560-64, 
1977) or enzymatic (F. Sanger et al., Proc. Natl. Acad. Sci. USA 
74:5463-67, 1977) methods. 

For the resolution and detection of target nucleic acids any mass 
5 determination method, such as, but are not limited to, chromatography, 
gel electrophoresis, capillary zone electrophoresis and mass spectrometry, 
is used. Mass spectrometric formats, include, but are not limited to, are 
matrix assisted laser desorption ionization (MALDI), electrospray (ES), ion 
cyclotron resonance (ICR) and Fourier Transform. For ES, the samples, 

10 dissolved in water or in a volatile buffer, are injected either continuously 
or discontinuously into an atmospheric pressure ionization interface (API) 
and then mass analyzed by a quadrupole. The generation of multiple ion 
peaks which can be obtained using ES mass spectrometry can increase 
the accuracy of the mass determination. Even more detailed information 

15 on the specific structure can be obtained using an MS/MS quadrupole 
configuration 

In MALDI mass spectrometry, various mass analyzers can be used, 
e.g., magnetic sector/magnetic deflection instruments in single or triple 
quadrupole mode (MS/MS), Fourier transform and time-of-flight (TOF) 
20 configurations as is known in the art of mass spectrometry. For the 

desorption/ionization process, numerous matrix/laser combinations can be 
used. Ion-trap and reflectron configurations can also be employed. 
Pair-matched nucleotide-based methods 
Forced Mass Modulation can be implemented using a 
25 deoxynucleotide set in which the mass of each base pair is identical, 

termed a pair-matched nucleotide set. A pair-matched nucleotide set can 
easily be formed, for example, by replacing dG (329.2 Da) in the set of 
naturally occurring nucleotides with 7-deaza-dG (328.2 Da). This forces 
the mass of each base pair to be 617.4 daltons: 
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(dA + dT) = (313.2 + 304.2) = 617.4 Da 
(dC + 7-deaza-dG) = (289.2 + 328.2) = 617.4 Da 
Many other pair-matched sets are possible using available nucleotide 
analogs. For this embodiment, the target DNA sequence can be 
5 composed entirely of the pair-matched nucleotide set. This can be 

accomplished by amplifying the target DNA sequence by PCR using the 
pair-matched nucleotide set prior to the sequencing reaction. 

A further requirement for this embodiment of Forced Mass 
Modulation is that the mass of each terminating base pair is unique and 
10 resolvable. The standard dideoxy terminators therefore cannot be used 
with the pair-matched nucleotide set described above, because the 
masses of all terminating base pairs are identical at 601 .4 daltons, except 
ddG:dC, which is 602.4 daltons. For the sake of clarity in this example, 
it is assumed that a set of mass-matched terminators is available ("ddM," 
15 defined as set of chain-terminating nucleotides that have exactly the same 
molecular weight ddA = ddC = ddG = ddT). If the mass of ddM is 
arbitrarily chosen to be 500 daltons, then the masses of the terminating 
base pairs are as follows: 

Terminatinp Base Pair Mass (Da) 

20 ddM: dC 789.2 

ddM: dT 804.2 
ddM: dA 813.2 
ddM: 7-deaza-dG 828.2 
In practice it is also possible to implement Forced Mass Modulation using 
25 a set of terminators that have different masses, this is discussed in detail 
below. 

Exemplary embodiments in which the mass shift is obtained using 
pair-matched nucleotides, where the mass of each nucleotide base-pair is 
selected so that the masses of all pairs are identical, are described in 
30 Example 4. In one embodiment thereof, the following steps are 
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performed: (i) the target nucleic acid is copied or amplified by a method 
such as PCR in the presence of the pair-matched nucleotide set prior to 
the sequencing or detection reaction; <ii) the target nucleic acid is 
denatured, and a partially duplex hairpin primer is annealed and ligated to 
5 the single-stranded template; (iii) the primer is extended in the presence of 
chain terminating nucleotides and pair-matched nucleotides to produce 
extension products; (iv) the masses of the extension products follow a 
periodic distribution that is determined by the mass of the pair-matched 
nucleotide set, and, (v) the target nucleic acid is detected by virtue of its 

10 molecular weight or its sequence is determined from the mass shift of 
each extension product from its corresponding periodic reference mass. 

In embodiments described above, the extending bases are pair 
matched, the mass of each terminating base pair is unique and 
resolvable, so that the mass shifts corresponding to each terminating 

15 base pair are unique. The nucleotide terminators are optionally mass- 
matched or can be of distinct masses as long as distinct values of mass 
shift are obtained for each terminating base pair. 

In another embodiment, the extension products are treated to 
produce blunt-ended double-stranded extension products by methods 

20 known to those of skill in the art, such as the use of single-strand specific 
nucleases. In an aspect of this embodiment, a plurality of target nucleic 
acids can be multiplexed in a single reaction by annealing each target 
nucleic acid to a primer of distinct molecular weight. The primers can be 
selected to differ in molecular weight by a value that is greater than the 

25 maximum mass shift, i.e., the difference in molecular weight between the 
heaviest and the lightest nucleotide terminating base pairs. Since double 
stranded nucleic acid can be analyzed, the effective sequence read is 
halved relative to the embodiment employing mass-matched nucleotides, 
but the number of molecules that can be multiplexed is doubled, due to 

30 the increase in period (the value of the mass of a base pair, rather than a 
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single mass-matched nucleotide). In exemplary embodiments, about 14 
to about 50 sequences are multiplexed. In detection embodiments, about 
50 or more molecules can be simultaneously multiplexed since only a 
single terminating base pair is added in the extension reaction. 

In another embodiment, the chain termination reactions can each 
be carried out separately using a standard nucleotide terminator, pair- 
matched nucleotides, and mass-labeled primers, if modified nucleotide 
terminators which are either mass-matched or provide distinct values of 
mass shift for each terminating base pair are not available. The reactions 
can be pooled prior to detection or sequence analysis. In one 
embodiment, the mass-labeled primers can have distinct values of 
molecular weight that give rise to unique values of mass shift or 
positional mass difference for each terminating base. 
Optimizing the mass spectrometry analysis of oligonucleotide mixtures 

In another method provided herein, nucleotide analogs are used to 
restrict the possible values of molecular weights that an oligonucleotide 
can possess relative to other oligonucleotides of the same length. The 
nucleotide analogs can be incorporated into the oligonucleotides using 
any suitable method, such as automated DNA synthesis (Oligonucleotides 
and Analogues: A Practical Approach, F. Eckstein, ed., IRL Press Oxford, 
1991) or by enzymatic replication using a polymerase and the requisite 
nucleotides and nucleotide analogs. 

For example, any two oligonucleotides with the same length n' 
with different base compositions can either 1 ) have exactly the same 
average molecular weight, or 2) have molecular weights no closer than a 
minimum value called the peak separation. In most cases, the peak 
separation will be a positive integer greater than one, but fractional values 
are theoretically possible. 

To illustrate an exemplary implementation of this method, the 
average molecular weight of the nucleotide analog 7-deaza-dG (328.2 
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daltons) can be substituted for Mg, into equation (ii) above, which defines 
Mavg as a function of the length "n" of the oligonucleotides in bases, as 
follows: 

(ii) Mavg = k + n'Mc + t(MT - Mc) + 3(Ma - Mc) + g(Mc - Mc), where 
5 Mc, Mt, Ma, Mg are the average molecular weights of each of the four 
nucleotide bases (cytosine, thymine, adenine, guanine); c, t, a, g 
represent the number of each base present in the oligonucleotide, the sum 
thereof, i.e., c + t + a + g = n', the total oligonucleotide length in 
bases; and the term k is a constant representing the mass of any other 
10 chemical groups on the molecule, such as terminal phosphates. 

Substituting the masses of the naturally occurring bases dC, dT and dA in 
DNA (to one-tenth dalton), and of 7-deaza-dG, 

Mc = 289.2 

Mt = 304.2 

15 Ma = 313.2 

Mg = 328.2 

and following simplification, the equation reduces to: 

Mavg = k + 289. 2n' + 1 5t + 24a + 39g 
Extracting the common factor from the last three terms yields 
20 (vi) Mavg = k + 289. 2n' + (5t + 8a + 13g) x 3 

In this example, the minimum peak separation is three daltons. It is 
not possible to identify or detect two oligonucleotides of the same length 
with different molecular weights that are closer than three daltons. 
Oligonucleotides with average masses closer than three daltons the 
25 oligonucleotides are detected if they are of different lengths. 

As a second example, Mt can be substituted with the molecular 
weight of a hypothetical nucleotide analog whose mass is 305.2 into 
equation (ii), yielding 

Mavg = k + 289. 2n' + 16t + 24a + 40g 
30 Extracting the common integer factor from the last three terms yields 
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(vii) Mavg = k + 289. 2n' + (2t + 3a + 5g) x 8 
for a minimum peak separation of eight daltons. Thus, appropriate 
selection of nucleotide analogs permits construction of nucleotide sets 
that provides sufficient peak separation for adequate resolution by mass, 
5 such as in a time-of-flight mass spectrometer. The trade-off for a greater 
peak separation is a greater number of base compositions that have 
exactly the same mass for a given oligonucleotide length. The maximum 
number of allowed oligonucleotide masses, L, for a given nucleotide set, 
is given by 

10 (Vii'l) L = n'dVlheavv - Mlight ) + 1, 

S 

where n' is the oligonucleotide length in bases, S is the peak separation, 

Mtight the mass of the lightest nucleotide in the set, Mheavy is the mass of 

the heaviest nucleotide in the set. The number of allowed oligonucleotide 

15 masses scales in direct proportion to the base length and inversely with 

the peak separation, but not all possible mass values will be represented 

for a given oligonucleotide length, particularly for small n. The average 

density of different base compositions per allowed mass value, D, can be 

obtained by dividing equation (v) by (viii) 

20 D = Ntotal . 

L 

which expands into 

(ix) D = S(n' + l)(n' + 2(n' -I- 3) 

6(n'(M G - Mc) + S) 

25 using a typical nucleotide set with G as the heaviest base and C as the 
lightest. The density function scales in direct proportion to the peak 
separation and as a quadratic function of the oligonucleotide length in 
bases. In practice, the average density of base compositions per allowed 
mass value predicated by equation (ix) will be somewhat lower than the 

30 actual density of base compositions per observed mass value, because 
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not all allowed masses will always be represented. The Examples 
describe implementation of the methods for sequencing. 
System and Software method for Force Mass Modulation 

Also provided are systems that automate the methods for 
5 determining a nucleotide sequence of a target nucleic acid or the 
detection methods provided herein using a computer programmed for 
identifying the sequence or target nucleic acid identity based upon the 
methods provided herein. The methods herein can be implemented, for 
example, by use of the following computer systems and using the 

10 following calculations, systems and methods. 

An exemplary automated testing system contains a nucleic acid 
workstation that includes an analytical instrument, such as a gel 
electrophoresis apparatus or a mass spectrometer or other instrument for 
determining the mass of a nucleic acid molecule in a sample, and a 

15 computer capable of communicating with the analytical instrument (see, 
e.g., copending U.S. application Serial Nos. 09/285,481, 09/663,968 and 
09/836,629; see, also International PCT application No. WO 00/60361 
for exemplary automated systems). In an exemplary embodiment the 
computer is an IBM compatible computer system that communicates with 

20 the instrument using a known communication standard such as a parallel 
or serial interface. 

For example, systems for analysis of nucleic acid samples are 
provided. The systems include a processing stations that performs a 
forced mass modulation chain extension reaction; a robotic system that 

25 transports the resulting products from the processing station to a mass 
measuring station, where the masses of the products of the reaction are 
determined; and a data analysis system, such as a computer programmed 
to identify nucleotides using forced mass modulation data, that processes 
the data from the mass measuring station to identify a nucleotide or 

30 plurality thereof in a sample or plurality thereof. The system can also 



-40- 



25491 -2408B 



include a control system that determines when processing at each station 
is complete and, in response, moves the sample to the next test station, 
and continuously processes samples one after another until the control 
system receives a stop instruction. 
5 The computer can be part of the instrument or another system 

component or it can be at a remote location. A computer system located 
at a site distant from the instrument can communicate wit the instrument, 
for example, through a wide area network or local area communication 
network or other suitable communication network. The system with the 

10 computer is programmed to automatically carry out steps of the methods 
herein and the requisite calculations. For embodiments that use 
mass-matched deoxyriboucleotides, a user enters the primer sequence or 
primer mass, the periodic reference mass and mass of an individual 
mass-matched deoxyonucleotide. These data can be directly entered by 

15 the user from a keyboard or from other computers or computer systems 
linked by network connection, or on removable storage medium such as a 
CD-ROM, minidisk (MD), DVD, floppy disk or other suitable storage 
medium. Next the user causes execution software that operates the 
system in which the mass spectrum of the extension products is 

20 generated. The Forced Mass Modulation software performs the steps of 
obtaining the masses of the fragments generated by the sequencing 
reaction and measured by the analytical instrument, and determining the 
identity of a nucleotide at any base position or the positional mass 
difference. The identity of the nucleotide at each base position is 

25 determined by comparing the calculated Mdiff[n] values to a database of 
previously calculated values of Mdiff for each of the chain terminating 
nucleotides. 

Mdiff[n] = Mobs[n] - MpR[n], 

where: 

30 (i) MpR[n] = (Mprimer + Mlight) + (n - 1 ) Phase, 
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in which n is the base position, MpR[n] is the n th periodic reference mass, 
Mprimer is the mass of the primer, M light is the mass of the lightest 
nucleotide terminator and Pbase is the base periodicity in daltons. The 
observed masses of the sequencing reaction products are given by the 
5 following equation: 

(ll) Mobs[n] = Mprimer + (n - 1 ) Pbase + Mterm[n], 

where n is the base position, Mobsfn] is the n th observed mass, Pbase is the 
base periodicity, and Mtermln] is the mass of the n th terminating nucleotide 
in daltons. The positional mass differences for the sequence can be 
10 obtained by subtracting equation (i) from equation (ii) and evaluating at 
every base position n: 

where Mdiff[n] is the n th positional mass difference. This relation simplifies 
to: 

(iii) Mdiff[n] = Mterm[n] - Mlight. 

15 Hence, the periodicity is determined by the mass of the mass- 

matched nucleotide and the shift is the difference in location of a peak 
resulting from the chain terminator. For example, in Figure 2, the lightest 
terminator is ddC, and the differential is O for C, 40 for G, 34 for A, 1 5 
for T. The selected mass matched nucleotide has a mass of 310 Da. The 

20 primer in Figure 2a has a mass of 3327 Da and the first peak would be at 
3600 if the first nucleotide in the extension product were C (0 shift). 
Since the first peak is at 3640, the shift is 40 Da. Therefore the first 
nucleotide is G, corresponding to a shift from the periodicity of 310 Da 
generated by the mass-matched nucleotides. 

25 Detection methods 

The methods herein may be used with any method for detection of 
nucleic acids based on molecular mass known to those of skill in the art, 
particularly methods in which a primer is extended. Such methods are 
modified by extending using mass matched nucleotides and/or chain 

30 terminators in extension reactions. Alternatively, or additionally, 
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amplification reactions may be performed using mass-matched nucleotides 
or pair-matched sets of nucleotides. These methods can be readily 
multiplexed using the methods and nucleic acid molecules provided 
herein. 

5 Detection methods and protocols, including those that rely on mass 

spectrometry (see, e.g., U.S. Patent No. 6,194,144; 6,225,450; 
5,691,141; 5,547,835; 6,238,871; 5,605,798; 6,043,031; 6,197,498; 
6,235,478; 6,221,601; 6,221,605; International PCT application No. WO 
99/31273, International PCT application No. WO 98/20019), can be 
10 modified for use with the methods herein by using mass-matched 
nucleotides for extension or pair matched duplexes for hybridization 
reactions. 

Among the methods of analysis herein are those involving the 
primer oligo base extension (PROBE) reaction with mass spectrometry for 

15 detection. In such reactions, the primer will be extended by mass- 
matched nucleotides. The methods herein are designed for multiplexing 
so that a plurality of different primers can be extended at different loci in 
the same reaction. The PROBE method uses a single detection primer 
followed by an oligonucleotide extension step to give products, which can 

20 be readily resolved by mass spectrometry, and, in particular, MALDI-TOF 
mass spectrometry. The products differ in length depending on the 
presence or absence of a polymorphism. In this method, a detection 
primer anneals adjacent to the site of a variable nucleotide or sequence of 
nucleotides and the primer is extended using a DNA polymerase in the 

25 presence of one or more dideoxy NTPs and, optionally, one or more 
deoxy NTPs. The resulting products are resolved by MALDI-TOF mass 
spectrometry. The mass of the products as measured by MALDI-TOF 
mass spectrometry makes possible the determination of the nucleotide(s) 
present at the variable site. Use of primers containing mass-matched 
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bases increases the resolving power of the reaction and permit 
simultaneous detection of a plurality of mutations (polymorphisms). 

These methods can be automated (see, e.g., copending U.S. 
application Serial No. 09/285,481 and published International PCT 
5 application No. PCT/US00/081 1 1, which describes an automated process 
line) and performed in a system that includes a computer programmed for 
analysis of the mass data as described above. 

The analyses can be performed on chip based formats in which the 
target nucleic acids or primers are linked to a solid support, such as a 

10 silicon or silicon-coated substrate, preferably in the form of an array. 
Generally, when analyses are performed using mass spectrometry, 
particularly MALDI, small nanoliter volumes of sample are loaded on, such 
that the resulting spot is about, or smaller than, the size of the laser spot. 
It has been found that when this is achieved, the results from the mass 

15 spectrometric analysis are quantitative. The area under the signals in the 
resulting mass spectra are proportional to concentration (when normalized 
and corrected for background). Methods for preparing and using such 
chips are described in U.S. Patent No. 6,024,925, co-pending U.S. 
application Serial Nos. 08/786,988, 09/364,774, 09/371,150 and 

20 09/297,575; see, also U.S. application Serial No. PCT/US97/201 95, 

which published as WO 98/20020. Chips and kits for performing these 
analyses are commercially available from SEQUENOM under the 
trademark MassARRAY. MassArray relies on the fidelity of the enzymatic 
primer extension reactions combined with the miniaturized array and 

25 MALDI-TOF (Matrix-Assisted Laser Desorption lonization-Time of Flight) 
mass spectrometry to deliver results rapidly. It accurately distinguishes 
single base changes in the size of DNA fragments associated with genetic 
variants without tags. 
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The following Examples are included for illustrative purposes only 
and are not intended to limit the scope of the invention. 

EXAMPLE 1 

Forced Mass Modulation using Mass-Matched Deoxynucleotides 

For this implementation, a set of nucleotide analogs for the four 
bases in DNA are selected (Amersham Pharmacia Biotech) such that each 
base has exactly the same molecular weight, termed a mass-matched 
deoxynucleotide set. This is achieved by judiciously choosing chemical 
modifiers of the existing bases or by the using a universal base analog 
such as deoxyinosine, which can form base pairs with more than one of 
the natural bases. For this example, the mass of each deoxynucleotide 
("dN") in the mass-matched set has the arbitrarily selected value of 310 
daltons, but any other value suffices. The sequencing reaction is 
performed as follows: 1 ) a primer is annealed to the target to be 
sequenced; 2) the resulting structure is subjected to an extension reaction 
using a suitable polymerase in the presence of the mass-matched 
nucleotide set and the four standard dideoxynucleotide terminators. The 
products and molecular masses of such a reaction are shown with a 
simulated mass spectrum in Figure 2a. The base periodicity is the mass 
of dN, or 310 daltons. The identity of a nucleotide at any base position is 
given by the positional mass difference, defined as the distance in daltons 
between the observed peak and the nearest periodic reference mass, 
which occurs every 310 daltons. In this example, the first periodic 
reference mass is defined as the (primer mass + ddC), or (3327 + 273) 
= 3600 daltons. The second periodic reference mass would be 3600 
plus the base periodicity or (3600 + 310) = 3910, and so on. 
Expressed in terms of the base position n: 

(i) MPR[n] = (Mprimer + Mlight) + (n - 1 ) Phase, 

where n is the base position, MpR[n] is the n th periodic reference mass, 
Mprimer is the mass of the primer, Miight is the mass of the lightest 
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nucleotide terminator and Pbase is the base periodicity in daltons. The 
observed masses of the sequencing reaction products are given by the 
following equation: 

(ii) Mobsfn] = Mprimer + (n - 1 ) Pbase + Mterm[n], 

5 where n is the base position, Mobsfn] is the n th observed mass, Pbase is the 
base periodicity, and Mterm[n] is the mass of the n th terminating nucleotide 
in daltons. The positional mass differences for the sequence can be 
obtained by subtracting equation (i) from equation (ii) and evaluating at 
every base position n: 
10 Mdiff[n] = Mobs[n] - MpR[n], 

where Mdiff[n] is the n ,h positional mass difference. This relation simplifies 
to: 

(lii) Mdiff[n] = Mterm[n] - Mlight. 

Inspection of equation (iii) reveals that Mdiff can only take on four distinct 
15 values, each corresponding to a different nucleotide terminator: 
Mdiff ["ddC"] = (273.2 - 273.2) = 0 
Mdiff ["ddT"] = (288.2 - 273.2) = 15 
Mdiff["ddA"] = (297.2 - 273.2) = 24 
Mdiff ["ddG"] = (313.2 - 273.2) = 40. 
20 Hence, the identity of the nucleotide at every base position in the target 
sequence can be determined by comparing each calculated positional 
mass difference with the values in the table above. Since the values that 
Mdiff can assume depend only on the choice of nucleotide terminators 
used in the sequencing reaction, it is possible to tailor the positional mass 
25 differences so that they are resolvable for any particular mass 

spectrometer. For example, replacing the terminator ddT with its analog 
5-bromo-dideoxyuridine (353.1 daltons) yields a positional mass 
difference of (353.1 - 273.2) = 79.9 Da for termination at T positions in 
the target sequence. This type of nucleotide substitution can be 
30 particularly valuable for lower-resolution mass spectrometers, as it 
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possible to maintain the sequence read accuracy without requiring any 
additional mass spectra. 

Further inspection of equation (iii) reveals that each observed mass 
value can be at most 40 daltons heavier than the nearest periodic 
5 reference mass. This limit is termed the maximum mass shift and is 
defined as the mass difference between the heaviest nucleotide 
terminator and the lightest. Resolving a second target sequence by 
Forced Mass Modulation with the standard dideoxy terminators is possible 
in a single spectrum so long as the primer for the second sequence is at 

10 least 40 daltons heavier (the maximum mass shift) than the primer for the 
first sequence, thus insuring that the peaks for each sequence never 
overlap in mass. 

In practice, it is recommended most mass spectrometric formats 
that the second primer is at least about 60 daltons heavier than the first 

15 primer, as each observed peak will have a finite width. Figure 2b shows 
a target second sequence resolved on the same mass spectrum shown in 
Figure 2a, using a primer heavier by 77 daltons. The peaks corresponding 
to the reaction products from the first target sequence can fall within the 
shaded regions of the spectrum in Figure 2b, which can never intersect 

20 peaks from the second target sequence. Unambiguous resolution of both 
sequences is possible in this arrangement because each peak can be 
uniquely assigned to a nucleotide, a base position, and a target sequence. 
This method is designated Mass Spectrum Division Multiplexing herein, 
and it is implemented using mass-staggered primers. Figure 3 shows four 

25 different sequences resolved in a single spectrum using a set of 

mass-staggered primers that are separated in mass by integer multiples of 
77 daltons (77, 154, and 231 daltons). 

The theoretical upper limit on the number of sequences that can be 
multiplexed in a single mass spectrum is given by the following equation: 

30 
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(iv) L' = Pbase 
Smax, 

where L' is the upper limit, Pbase is the base periodicity, and Smax is the 
maximum mass shift in daltons. For the nucleotide set and terminators 
5 used in this example, L = (310 / 40) = 7.75, or approximately seven. 
Increasing the number of sequences that can be multiplexed in a single 
spectrum, can be achieved by implementing one or both of an increase in 
the base periodicity, and a reduction of the maximum mass shift. The 
base periodicity can be increased by choosing a mass-matched nucleotide 

10 set that has a higher molecular weight for dN. It is simpler to lower the 
maximum mass shift by careful use of the nucleotide terminators and their 
analogs. For example, if the sequencing reactions were performed using 
only the terminators ddC, ddT, and ddA, then the maximum mass shift 
becomes (mass of ddA - mass of ddC) = (297 - 273) = 24 Da. In this 

15 case the upper limit on the number of sequences that can be multiplexed 
is L = (310 / 24) = 12.92, or approximately twelve. In situations where 
complete sequence information is not required, such as diagnostic 
sequencing, a great reduction in the number of required spectra can be 
realized by using fewer than four nucleotide terminators. If the 

20 sequencing reaction is performed using only a single nucleotide 

terminator, the maximum mass shift becomes identically zero, and the 
number of sequences that can be multiplexed in a single spectrum is 
limited only by the absolute resolution of the mass spectrometer in 
question. If a given mass spectrometer has an absolute resolution of 1 2 

25 Da in the mass range of the sequencing reaction products, then the 

maximum number of sequences that can be multiplexed is given by L = 
(310/12) = 25.83, or approximately twenty-five. 
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EXAMPLE 2 

Forced Mass Modulation using Pair-Matched Deoxynucleotides 

Implementation of Forced Mass Modulation using pair-matched 
nucleotides is shown in Figure 4. The basic requirement for this method 
5 is that the sequencing reaction products can be analyzed as 

double-stranded structures. Briefly, the steps in the reaction are as 
follows: 1) A partially duplex hairpin primer with a 3' overhang and a 5' 
phosphate group is annealed and ligated to the single stranded target 
sequence. 2) The resulting partially duplex structure is subjected to a 

10 sequencing reaction using the pair-matched nucleotide set described 
above along with the set of mass-matched terminators (ddM). 3) The 
products from the sequencing reaction are exposed to a strict single 
strand-specific nuclease that results in the production of blunt-ended 
hairpin structures ready for analysis by mass spectrometry. Figure 5 

15 shows the products and molecular masses of the nuclease digestion along 
with a simulated mass spectrum. 

Because the reaction products are double-stranded, they are forced 
to assume a quasi-periodic distribution with a base periodicity of 617.4 
daltons. The shaded regions on the spectrum shown in Figure 5 indicate 

20 the allowed mass ranges that can be occupied by the reaction products. 
The first periodic reference mass is at 10360 Da, which is the mass of 
the fully duplex hairpin primer plus a ddM:dC base pair. Expressing the 
periodic reference masses in terms of the base position n yields: 

<X) MpR[n] = (Mduplex + Mlight + MddM) + (n - 1 ) X Pbase 

25 Where MpR[n] is the n th periodic reference mass, Mduplex is the mass of the 
fully duplex primer, Miight is the mass of the lightest deoxynucleotide in the 
target, Pbase is the base periodicity, and MddM is the mass of ddM in 
daltons. The observed masses of the sequencing reaction products are 
given by the following equation: 

30 (Xi) Mobsfn] = Mduplex + MddM + (n - 1) X Pbase + Mtarg[n], 
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where n is the base position, Mobs[n] is the n ,h observed mass, and Mtar g [n] 
is the mass of the n th nucleotide in the target sequence past the priming 
site in the 3' -> 5' direction. 

In contrast to the mass-matched nucleotide set implementation that 
5 provides the sequence complementary to the template strand read in the 
5' -> 3' direction, the pair-matched nucleotide set implementation 
described herein directly reads the template strand in the 3' -> 5' 
direction. The positional mass differences for this implementation are the 
same as those in Example 1 , except that the mass difference 

10 corresponding to a termination on dG is 39 as opposed to 40 daltons, 
because 7-deaza-dG is exactly one dalton lighter than dG. Since double 
stranded DNA can be analyzed for this method to work, the effective 
sequence read length is halved, although the number of sequences that 
can be multiplexed is doubled, due to the increase in the base periodicity. 

15 As a demonstration of Forced Mass Modulation implemented 

without using mass-matched terminators, the positional mass differences 
for the above example using the following set of nucleotide terminators is 
calculated as follows: 



Terminator 


Nucleotide Analog 


Mass 


Base 
Pairing 


Mass of 
Base Pair 


T 


5-Bromo-dideoxyuridine 


353.1 


5-Br-ddU:dA 


666.3 


C 


5-Methyl-dideoxycytidine 


287.2 


5-Me-ddC:7-deaza-dG 


615.4 


A 


Dideoxyadenosine 


297.2 


ddA: dT 


601.4 


G 


Dideoxyinosine 


298.2 


ddl: dC 


587.4 



25 The positional mass difference at every base position is given by: 
(xii) Mditf[n] = M P air[n] - Mnghtest, 
where Mdiff[n] is the n ,h positional mass difference, M P air[n] is the mass of 
n th terminating base pair, and Mnghtest is the mass of the lightest 
terminating base pair in daltons. Substituting in the values from the table 

30 above yields: 
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Mdiff["G"] = (587.4 - 587-4) = 0 
Mdifff'A"] = (601.4- 587-4) = 14 
Mdiff["C"] = (615.4- 587.4) = 28 
Mdiff["T"] = (666.3 - 587-4) = 78.9 
5 Since each terminating base pair has a unique positional mass difference, 
the base sequence can be determined unambiguously. The maximum 
mass shift in this case is 78.9 daltons. When choosing a set of 
terminating nucleotides it is important to select the set such that the 
positional mass difference for each base termination is distinct and 
10 resolvable by mass. 

If modified nucleotide terminators are not used, it is still possible 
to implement Forced Mass Modulation by carrying out each of the four 
termination reactions separately using mass-labeled primers rather than 
modified terminators, combining all reaction products, and then obtaining 
15 a mass spectrum. In order to produce the same positional mass 
differences as shown in Example 1 , using a set of pair-matched 
nucleotides and the standard dideoxy terminators, the following primer 
mass shifts are required: 

Termination Reaction Primer Mass 
20 C "reference " primer 

T reference primer + 15 Da 

A reference primer + 24 Da 

G reference primer + 39 Da 

This method is essentially equivalent to multiplexing four single-nucleotide 
25 sequencing reactions in the same spectrum, except that all the 

sequencing products originate from the same priming site but terminate 
on different nucleotides. 
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EXAMPLE 3 

Forced Mass Modulation in the Detection and Scoring of Single Nucleotide 
Polymorphisms 

Forced Mass Modulation can be used to simplify the analysis of 
5 closely related sequence variants, as is required in the detection and 
scoring of single nucleotide polymorphisms. Figure 6 shows three 
sequence variants that differ from each other only at a single base 
position sequenced by a conventional Sanger reaction. The mass 
distribution of the reaction products is so complex that it can be 
10 uninterpretable, even if the base sequences of the variants are known a 
priori. 

Figure 7 shows the same three variants sequenced by Forced Mass 
Modulation using mass-matched deoxynucleotides <dN = 310 Da) and the 
standard dideoxy terminators. The positions and identities of the 
15 single-nucleotide changes are immediately apparent from the mass 
spectrum. Since the masses of the sequencing reaction products are 
constrained to fall within the shaded regions of the spectrum in Figure 7b, 
it is possible to multiplex other sequences on the same spectrum. 
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EXAMPLE 4 

Base Composition Density Distributions for the Total Set of possible 7- 
base Oligonucleotides 

For this implementation, three sets of 7-base oligonucleotides 

5 comprising all possible base compositions for a 7-base oligonucleotide can 

be obtained; the first set comprising the four natural bases (dA, dG, dC 

and dT), the second set comprising three of the natural bases (dA, dC 

and dT) and the nucleotide analog 7-deaza-deoxyguanosine (7-deaza-dG) 

substituted for dG, and the third set comprising three of the natural bases 

10 (dA, dG and dC) and the nucleotide analog deutero-deoxythymine 
(deutero-dT) substituted for dT. Figure 8 shows the actual base 
composition density distributions for the total set of possible 7-base 
oligonucleotides using the three different nucleotide sets. Note that for 
the set of naturally occurring bases (Figure 8a), nearly every base 

15 composition has its own distinct mass value, but most of these mass 
values are spaced only one dalton from each other. Increasing the peak 
separation to three daltons by substitution of dG with 7-deaza-dG (Figure 
8b) markedly increases the average number of base compositions per 
observed mass, particularly for those masses in the center of the range, 

20 but any two oligonucleotides of the same length with different molecular 
weights will have to be separated by at least three daltons. Similarly, 
substitution of dT with deutero-dT (Figure 8c) gives a minimum peak 
separation between oligonucleotides having the same length but different 
molecular weights of eight daltons. The trade-off for a greater peak 

25 separation is a greater number of oligonucleotides that have exactly the 
same mass for a given oligonucleotide length. 



Since modifications will be apparent to those of skill in this art, it is 
intended that this invention be limited only by the scope of the appended 
30 claims. 
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